Neutrino production in UHECR proton interactions in the infrared 
background 

Todor Stanev 

Bartol Research Institute, University of Delaware 
Newark, DE 19716 USA 

We discuss the contribution of proton photoproduction interactions on the isotropic infrared/optical background 
to the cosmic neutrino fluxes. This contribution has a strong dependence on the proton injection energy spectrum, 
and is essential at high redshifts. It is thus closely correlated with the cosmological evolution of the ultra high 
energy proton sources and of the inrared background itself. These interactions may also contribute to the source 
fluxes of neutrinos if the proton sources are located in regious of high infrared emission and magnetic fields. 



The assumption that the Ultra High Energy 
Cosmic Rays (UHECR) are nuclei (presumably 
protons) accelerated in luminous extragalactic 
sources provides a natural connection between 
these particles and ultra high energy neutrinos. 
This was first realized by Berezinsky&Zatsepin pQ 
soon after the introduction of the GZK effect [2] ■ 
The first realistic calculation of the generated 
neutrino flux was made by Stecker |3] . The prob- 
lem has been revisited many times after the paper 
of Hill&Schramm 0] who used the non-detection 
of such neutrinos to limit the cosmological evolu- 
tion of the sources of UHECR. 

These so called cosmological neutrinos are 
produced in photoproduction interactions of 
the UHECR with the ambient photon fields, 
mostly with the microwave background radiation 
(MBR). The GZK effect is the limit on the highest 
energy a cosmic ray proton can retain in propa- 
gation through the MBR. It sets a cutoff in the 
cosmic ray energy spectrum in case the UHECR 
sources are isotropically and homogeneously dis- 
tributed in the Universe. The physics of these 
photoproduction interactions is very well known. 
Although the energy of the interacting protons 
is very high, the center of mass energy is low, 
mostly at the photoproduction threshold. The in- 
teraction cross section is studied at accelerators 
and is very well known. Most of the interactions 
happen at the A + resonance where the cross sec- 
tion reaches 500/ib. The mean free path reaches 



a minimum of 3.4 megaparsecs (Mpc) at proton 
energy of 6xl0 20 eV. The average energy loss of 
10 20 protons is about 20% per interaction and 
slowly increases with the proton (and center of 
mass) energy. 

The fluxes of cosmological neutrinos are, how- 
ever, very uncertain because of the lack of cer- 
tainty in the astrophysical input. The main pa- 
rameters that define the magnitude and the spec- 
tral shape of the cosmological neutrino fluxes 
are: the total UHECR source luminosity Lqr, 
the shape of the UHECR injection spectrum 
olcr in the case of power law spectrum, the 
maximum UHECR energy at acceleration E max 
and the cosmological evolution of the UHECR 
sources. These are the same parameters that 
Waxman&Bahcall J5] used to set a limit on the 
neutrino fluxes generated in optically thin sources 
of UHECR. 

The microwave background is not the only uni- 
versal photon field that has to be taken in con- 
sideration. Especially interesting is the isotropic 
infrared and optical background (IRB). The num- 
ber density of IRB is smaller than that of MBR 
by more that two orders of magnitude. On the 
other hand, protons of lower energy can interact 
on the IRB, and the smaller number density has 
to be weighted with the larger flux of interacting 
protons. The present Universe is optically thin 
to 10 19 eV and lower energy protons, but even at 
small redshift the proton interaction rate quickly 
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increases. This is different from the interactions 
on MBR, where the interacting protons quickly 
lose their energy even at z=0. The cosmological 
evolution of UHECR injection is thus of major 
importance for the contribution of such interac- 
tions to the flux of cosmological neutrinos. 

We use the IRB model of Franceschini et al [H] 
shown in Fig. together with the MBR in terms 
of energy density. The model consists of two com- 
ponents: 'star', near infrared, which covers the 
higher photon energies, and 'dust', far infrared 
that continues down to MBR. The total IRB num- 
ber density is significantly smaller than that of 
MBR. The model yields 1.6 photons/cm 3 , a factor 
of 250 less than the MBR. The IRB is measured 
directly after subtraction of point sources and is 
also estimated from the absorption of TeV pho- 
tons coming from extragalactic sources ■ These 
estimates affect mostly the near infrared part of 
the spectrum. Photons of wavelength above 40 
fim affect only the 7-ray fluxes above 10 TeV 
where the statistics is usually low and the flux 
decrease could also be due to absorption in the 
7-ray sources. 
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Figure 1. The energy density in MBR and IRB 
according to the model of Ref. UJ. The data 
points are from analyses of the DIRBE measure- 
ments piTU] . 



In addition to the lower total photon density 



the IRB covers much wider wavelength range 
than the microwave background, and its photon 
density per unit energy is even smaller. The in- 
teractions of UHECR on IRB photons are indeed 
very rare in the present universe. Fig.[5]shows the 
fraction of the proton energy that is converted to 
neutrinos as a function of the proton energy in 
propagation on a distance of 200 Mpc. 
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Figure 2. Fraction of the proton energy that 
is converted to neutrinos on propagation on 200 
Mpc. The histogram shows interaction on MBR 
only, while the points represent interactions of the 
photon spectrum shown in Fig. 



In the derivation of the neutrino limit Wax- 
man&Bahcall use cosmic ray source luminosity 
Lqr = 4.5 ± 1.5 x 10 44 erg/Mpc 3 /yr between 
10 19 and 10 21 eV for power law cosmic ray en- 
ergy spectrum with a — 2. The assumption is 
that no cosmic rays are accelerated above 10 21 
eV. The cosmological evolution of the source lu- 
minosity is assumed to be (1 + z) 3 to z = 1.9 
then flat to z—2.7 with an exponential decay at 
larger redshifts. We will first use the parameters 
of this limit to find the contribution of the proton 
interactions on IRB. 

The resulting + spectrum for a cosmo- 
logical model with Qa = 0.7, £1m = 0.3 and Ho 
= 75 km/s/Mpc is shown with a dotted line in 
Fig. |3 The flux peaks at 10 16 3 eV at 2.5 xl0~ 18 
cm -2 s -1 ster _1 . The peak is at energy lower than 
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the peak of the MBR interactions (shown with a 
dash-dot line) by a factor of 20, and its magnitude 
is also lower by a factor of 10. Next we show in the 
same figure with a dashed line the contribution of 
IRB for a scenario in which the injection spectral 
index is changed to a = 2.5 and all other param- 
eters are the same. There is a noticeable shift 
of the peak position to still lower energy. The 
peak is now located at 10 15 7 eV and is higher by 
a factor of about 7. The contribution of IRB is 
now smaller than that of MBR (a — 2) only by 
about 30%. The highest curve in Fig. shows 
the IRB contribution for a = 2.5 and cosmolog- 
ical evolution with n = 4 and then constant to 
z=10 followed by an exponential decrease. The 
location of the peak does not change but its mag- 
nitude increases by almost a factor of three. It is 
now 50% higher than the 'standard' MBR gener- 
ated cosmological neutrinos. It is obviously not 
correct to compare fluxes obtained with different 
assumptions for the cosmological evolution and 
we do it only to have a feeling for the magnitude 
of the neutrino fluxes. The a = 2.5 spectra de- 
crease the flux of cosmological neutrinos of energy 
above 10 19 eV. 

Both the spectral shape and the cosmological 
evolution of the UHECR sources affect the con- 
tribution of the IRB to the cosmological neutrino 
flux. The most important factor, however, is the 
shape of the injection spectrum. It is worth to 
note that the maximum proton energy at accel- 
eration does not affect the IRB generated fluxes, 
since they are due mostly to protons of energy 
below 10 20 eV, as can be observed in Fig. [3 

At energy about 3xl0 18 eV the cosmological 
fluxes of + are very close to the limit for 
source neutrinos. The reason is simple - in prop- 
agation from large distances protons lose almost 
all of their energy in interactions on MBR. An 
interesting feature is the flux of v e (not shown), 
which peaks at energy about 3xl0 15 eV. The ori- 
gin of this flux is neutron decay, and a small v e 
flux is generated in neutron interactions on MBR. 

The cosmological evolution of the sources 
(n=3) increases the fluxes by about a factor of 
five compared to a no-evolution scenario. The in- 
crease, however, is not energy independent |12) . 
The highest energy neutrinos are generated at 
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Figure 3. Fluxes of cosmological neutrinos 
[y^ + Dfj,) generated only by interactions on IRB. 
All three calculations use the UHECR luminosity 
derived by Waxman . The power law spectral 
indices and the cosmological evolution of UHECR 
sources n are given by each curve. The dash- 
dotted line shows the 'standard' n—3 cosmologi- 
cal neutrino flux from interactions in the MBR. 



small redshifts. The low energy neutrinos come 
from high redshifts because of two reasons: the 
threshold energy of protons for photoproduction 
interaction decreases, and the generated neutri- 
nos are further redshifted to the current epoch. 
The standard flux (a=2.0, n=3) would gener- 
ate about 0.4 neutrino induced showers per km 3 
year in the IceCube ^3] neutrino detector and 0.9 
events with energy above 10 19 eV in the Auger [T5| 
observatory (for target mass of 30 km 3 of water) 
assuming that at arrival at Earth the flavor ra- 
tio v e : : v T is 1:1:1 because of neutrino os- 
cillations. It is difficult to estimate the rate in 
EUSO because of its yet unknown energy 
threshold. These events come from the NC inter- 
actions of all neutrinos, CC interactions of v e , the 
hadronic (y) part of the CC interactions of muon 
and tau neutrinos and from r decay. Although 
very prominent, the Glashow resonance does not 
produce high rate of events because of its narrow 
width. Ice Cube should also detect very energetic 
muons with a comparable rate which is difficult 
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to predict without detector Monte Carlo simula- 
tions. 
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Figure 4. Comparison of the cosmological neu- 
trino fluxes with the Waxman&Bahcall limit, 
which is given as an shaded area for the 'stan- 
dard' power law injection spectrum and cosmo- 
logical evolution. The thick white line shows the 
limit derived in Ref. ^21- The dashed line shows 
the flux of cosmological neutrinos generated in 
interactions on MBR for the 'standard' parame- 
ters, and the solid one - for a — 2.5 and n = 4. 
The squares show the fluxes generated on the to- 
tal photon background shown in Fig.^ the open 
squares are for the 'standard' parameters and the 
full ones - for a = 2.5 and n = 4. 



Changing the proton injection spectrum to a 
power law with a = 2.5 moves the maximum 
of the cosmological neutrino flux to lower en- 
ergy and increases the contribution of the inter- 
actions on IRB. At the same time the flux of 
higher energy cosmological neutrinos decreases. 
The shower event rates in IceCube and Auger be- 
come 0.44 and 0.31 respectively. 

Assuming a stronger source evolution, (1 + z) 4 
makes a big difference in the expected fluxes. 
With a power law source spectrum with a = 2.5 it 
generates 1.2 events in IceCube and 0.66 events 
in Auger. The cosmological neutrino spectrum 
for + is shown with full squares in Fig. 0] 



The contribution of interactions on MBR is shown 
with a solid line. 

The biggest uncertainty in these results, which 
is not listed above, is the cosmological evolu- 
tion of the infrared/optical background. The es- 
timates above assume that it is the same as of 
MBR, i.e. that the IRB was fully developed at 
z=8, which is the limit of the redshift integration. 
This does not seem to be a realistic assumption, 
although models of the IRB emission ^7] predict 
very strong evolution of the far infrared emission, 
especially between redshifts of 10 to 100. 

The maximum proton energy at acceleration 
Emax is unknown, but having in mind the highest 
energy Fly's Eye shower of 3xl0 20 eV one should 
expect that astrophysical sources accelerate pro- 
tons at least to 10 21 eV. The injection spectrum 
is also not very well determined since the result 
of proton propagation depends on the UHECR 
source distribution. Attempts to derive the in- 
jection spectrum in the case of isotropic homoge- 
neous source distribution end up with injection 
spectra not flatter than E~ 2A power law |18I19| . 
. The extreme case is developed by Berezinsky et 
al. |20j who derive an a=2.7 injection spectrum. 

The luminosity required for the explanation of 
the observed events above 10 19 eV grows with the 
spectral index, and in the case of Berezinsky et 
al. becomes 4.5 x 10 47 erg Mpc _3 yr _1 . Such steep 
spectrum would generate only a small event rate 
for neutrinos above 10 19 eV and would enhance 
the IRB contribution. 

Expressed in terms of (1 + z) n the cosmologi- 
cal evolution of different astrophysical objects is 
observed to be between n = 3 and 4. A strong 
evolution with n = 4, as used above, may be too 
optimistic, but not entirely out of range. As seen 
from Fig. 0] strong cosmological evolution does 
not only increase the total flux, but moves the 
peak of the cosmological neutrino spectrum to 
somewhat lower energy. 

Finally, the cosmic ray source luminosity, which 
was normalized to the flux of UHECR at 10 19 
by Waxman could easily be higher or lower 
by half an order of magnitude. One can then 
assume a pessimistic IceCube shower event rate 
of 0.1 event per km 3 yr and an optimistic rate of 
4-5 events. 
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It is obvious that a detailed calculation of the 
flux of cosmological neutrinos should include the 
interactions on the infrared background. We plan 
to do that with a better model of the IRB cosmo- 
logical evolution and describe the calculation in 
more detail in a forthcoming paper. One should 
also keep in mind that if the UHECR sources 
are located in regions of high infrared and opti- 
cal photon density, the fluxes of source neutrinos 
could increase. The effect may be much stronger 
if 10 19 eV and lower energy protons are contained 
in the region by high magnetic fields. 
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